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Abstract We give a direct tensor decomposition for any density matrix into 
Hermitian operators. Based upon the decomposition we study when the mixed 
states are separable and generalize the separability indicators to multi-partite 
states and show that a density operator is separable if and only if the separable 
indicator is non-negative. We then derive two bounds for the separable indicator 
in terms of the spectrum of the factor operators in the tensor summands. 
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1 Introduction 
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In the last decade quantum entanglement has played a remarkable role in many applications and 
become one of the key resources in the rapidly expanding fields of quantum information and quan- 
tum computation, especially in quantum teleportation, quantum cryptography, quantum dense 
coding and parallel computation [1, 2, 3]. A quantum state or density matrix is separable (or 
not entangled) if it is a convex sum of tensor product of quantum states. In this case the sepa- 
rable quantum state can be prepared in several different locations. There are two aspects in the 
question regarding quantum entanglement: the first is to judge whether a general quantum state 
is entangled or not, and the second is to answer how much entanglement remained after some 
noisy quantum process. In the case of pure states, the Bell inequality provides a useful tool to 
(^ ', tell separability from entanglement [4]. In [5, 6, 7, 8] separability problem was examined and 

important criteria were proposed from several viewpoints for the far more difficult case of mixed 
states including the PPT criterion and the range equality condition. In terms of measurement of 
entanglement other methods have been found, e.g. formation of entanglement [9, 10] and purifi- 
cation of formation [11, 12]. Recently further important and interesting works [13, 14, 15] have 
also been devoted solely to quantum entanglement and some criteria were proposed accordingly, 
in particular, [16] gives an operational and geometric approach to pairwise entanglement of two 
and three-dimensional composite quantum systems. Despite these important developments the 
question of separability still remains unsolved and is notoriously famous for its difficulty. 

Among the approaches to quantum separability it is highly needed to have an operational 
method to decompose the quantum states as tensor product. Such an idea was first studied in 
[17], where some necessary constraints were found to ensure an optimal separable approximation 
to a given density matrix, and then a numerical method is proposed to locate the optimal separable 
state for two-partite mixed states. In [18] a new algebraic mechanism was introduced to study 
the separability question for two partite mixed states. The idea was first to decompose the mixed 
density matrix as a summation of tensor products of Hermitian operators, and then we rearrange 
the sum to get the indicator. It was proved that the density matrix is separable if and only if the 
indicator is non-negative. Thus the indicator provides a new measurement for the separability. 
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In this article we will generalize this method to multi partite density operators. We will give 
a new operational method to decompose the density matrix as a summation of tensor products of 
Hermitian operators. Our new method at the simplest case is the fundamental fact that any 4x4 
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Hermitian operator is a span of composite Pauli spin matrices &i <S> Pj, where ctq 



ai = {°i i)' CT2 = oJ' and(T3= (J -°i 

Once the decomposition into tensor product is known, the idea of separability indicator [18] is 
generalized into multi-partite states and we show that the mixed states are separable if and only 
if the separability indicator is non-negative. In general it is hard to compute the separability indi- 
cator. For this purpose we provide several bounds, and hope that they will help in determination 
of the separability. 



2 Basic notions 

Let ifi(resp. i^be an m(resp. n)-dimensional complex Hilbert space, with \i),i = l,2...,m 
(resp. \j),j = l,2...,n) as an orthonormal basis. A bipartite mixed state is said to be separable 
if the density matrix can be written as 

i 

where < pi < 1, ^2pi = 1, p\ and p 2 are density matrices on H\ and H2 respectively. It is 

i 

a challenge problem to find such a decomposition or proving that it does not exist for a generic 
mixed state [5, 6, 7, 8]. 

We first introduce some notations. For an m x m block matrix Z with each block Zij of size 
n x n, i,j = l,2...,m. The realigned matrix Z is defined by 

Z = [wec(Zii),...,uec(Z m i),...,wec(Zi m ),...,wec(Z mm )]*, (2) 

where for any m x n matrix T with entries tij, vec(T) is defined to be 

Vec(T) = [tll,...,t m l,tl2,—,tl n ,...,t mn ] ■ 

Let A = A +y/—lA I be a complex Hermitian matrix, where A R and A 1 are real and imaginary 
parts of A. Let a be the canonical map from A to a real matrix: 

( AR Al \ 
a:A^^_ Al AR y (3) 

where A R and A 1 are the real and imaginary parts of A respectively. 

Let Q s be an m 2 x VHSUL — I matrix. If we arrange the row indices of Q s as 

{11, 21, 31,..., ml, 12, 22, 32,. ..,m2, ...,mm}, 

then all the entries of Q s are zero except those at 21 and 12 (resp. 31 and 13,...) which are 1 and 
-1 respectively in the first (resp. second,...) column. In other words, 

Qs = [{^21, — ei2J; {e3i, — ei3J; ...; \e mm -\, — e m _i ;Tn }], 



where {e2i, — ei 2 } is first column of Q s , with 1 and —1 at the 21 and 12 rows respectively; while 
{e 3 i, — ei3J is second column of Q s , with 1 and —1 at the 31 and 13 rows respectively; and so on. 



Let Q a be an m? x 



m(m+l) 



matrix such that 



Qa = [{en}; {e 2 i, ei 2 }; {e 3 i, ei 3 }, ...; {e 22 }; {e 32 , e 23 }, {e42, e 24 }; •••; {e m , m _i, e m _i, m }, {e mm }], 

where {en} is the column vector with 1 at the row m and zero elsewhere, and {eij,e±j} is the 
column vector with 1 at the ijth and jith. rows and zero elsewhere. 



Qi can be expressed as 



Qi 



Q s 
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where Q s and Q a are obtained by normalizing each column of Q s and Q a 
By replacing the dimension m with n, we have Q 2 . 
As an example we have for m=2 
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3 The tensor product decompositions of Hermitian matrices 



Let A be a Hermitian matrix on Hilbert space Hi <S> H 2 . In [18] we gave an operational method 
to decompose A as a tensor product of Hermitian matrices on Hi and i? 2 . We will give another 
method to decompose A and then generalize to the case of multi-tensor products. 

Let's recall the decomposition method in [18]. We express the matrix A in terms of real and 
complex parts: A = A + iA and realign both A and A into A and A respectively as in Eq. 
(2). Then we write 

' A* Ai\ _( An A\ 2 
-A 1 A* j ^ 2 I Mi A 22 



Q\ 



(4) 



Proposition 1. Let A be an mn x mn Hermitian matrix as rewritten in Eq. (4). Suppose the 
singular value decomposition of ^4 22 is ^4 22 = Ya=i VKuivj, where r is the rank and Aj(i=l,2...,r) 
are the non-zero eigenvalues of ^4 22 A 22 , and U{ (resp. vi) are the eigenvectors of the matrix 

A 22 A 22 (resp.A 22 A 22 ). Set Bi = y/Xiui, Ci = —Vi. Then we can decompose A as a tensor 
product 

r 



i=l 



where the m x m Hermitian matrices Bi 
ci + \J —\Ci are given by 



hi + yJ—V&i and the n x n Hermitian matrices Ci = 



vec(bi) 
-vec{Bi 



Qi 



o 

-Bi 



vec(ci) 
vec(Ci) 



Q-2 




Ci 



(5) 



The above result gives a constructive or operative method to decompose A as a tensor product. 
The existence of tensor decomposition has a simpler explanation. In fact, we know that the set of 
n x n Hermitian matrices is a real vector space of dimension n 2 , thus the dimension of Hermitian 
matrices of size ran x mn is exactly equal to the product of the dimension of size m x m and 
that of size m x m, hence the subspace of tensor product of Hermitian matrices of size n x n and 
that of size m x m must equal to the space of all Hermitian matrices of size mn x mn, which 
guarantees the existence. 

We observe that in general the space of real symmetric (antisymmetric) matrices can not be 
decomposed into tensor product of symmetric (antisymmetric) matrices. In fact, the difference 
between dimensions of the space of mn x mn symmetric matrices and that of the tensor product 
of symmetric matrices of size m x m and size n x n is 



mn + 1 
2 



m + 1 
2 



n+ 1 
2 



m 
2 



Similarly the difference between the dimensions of antisymmetric operators over C m x C r ' 
that of the tensor product of antisymmetric operators is 



and 



mn — 1 
2 



771—1 
2 



n — 1 
2 



m + 1 
2 



n+ 1 
2 



1. 



We can use induction to generalize Proposition 1 to multi-partite case. 

Theorem 1. Let A be an Hermitian matrix on space H\ (g H^ <g H% (g ... (g H n . A has tensor 
production decomposition like A = Yll=i B} <g -B 2 (g ... (g J3™, where £*, B|, ..., 1?" are Hermitian 
matrices on H±, H2, ■■■, H n respectively. 

We now present a practical method to decompose Hermitian matrices into tensor product of 
Hermitian matrices, thus giving a new constructive proof for Theorem 1. Let E^- be the unit 
square matrices of size n x n. If it is clear form the context, we will omit the superscript. To 
decompose the unit matrix E^ n , we write its indices i,j uniquely as follows: 
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i = (k — \)n + i', j - 
where 1 < k, I < m and 1 < i',f < n. Then we have 



(l-l)n + f, 



(6) 






E 



ki 



jpn 



(7) 



Equivalently we can picture the above decomposition as follows. We first view E^ n as an 
m x m block matrix with each entry as an n x n matrix. The resulted block matrix is still a 
unit-like matrix where all entries are zero except (k, Z)-entry, which is an n x n unity matrix itself, 
say Ei/ji. Then we immediately have EfJ n = E™ t (g) E™,-,. 

Example 1. Let (1 + 7b)~ l pi, be the density operator on C 2 <8) C 4 as follows. 
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We can decompose p\, by the above scheme. 



Pb 



H E n + E W + ^22 + -^27 + -^33 + -^38 + E M + E 61 + E 66 + E 72 + E 77 + -^83 



+ 



+ 



+ 



-^—(£55 + E 88 ) 



VT^W 



+ - — 2 ^ 58 + ^ 85 -' 

&(£?! ® 4 + £? 2 ® #12 + #11 ® #22 + #12 

^21 + -^22 



-^23 + -^11 ® ^33 + -^12 



-^34 



1 + 6 



-^22 + -^21 ® -^32 + -^22 



E^ + £^21 



v 33 



^3) 



(E, 



22 



# 1 + E. 



22 



i&) + 



VT^W 



(E? 



22 



£? 4 + EL ® E\ 



22 



^4i; 



For a different decomposition using the singular value decomposition, the reader is referred to 
[18]. 

This decomposition method can be generalized to Hermitian operators. Let A be a Hermitian 
matrix, then one can decompose A into a sum of real and imaginary parts: A = B + \J — 1C , 
where B (or C) is a symmetric (or antisymmetric) matrix. Let {Ey + Eji} be the basis for the 
symmetric matrices, and {Eij — Eji} be the basis for the antisymmetric matrices. It is enough 
to decompose the basis elements as tensor products of Hermitian matrices. Roughly speaking, 
one writes each basis element E^ ± Eji of size mn x ran as a block matrix, then transform it 
into a tensor product according to the position where the 1 or —1 appears. The main point is 
that we have to consider all Hermitian matrices to factor the basis elements (cf. the remark after 
Proposition 1). 

Specifically, by modulo n we write the indices i,j uniquely as in Eq. (6): i = i'(modn),j = 
f(modn) and k = [(i — l)/n] + 1,1 = [(j — l)/n] + 1. Here the representatives for Z n are taken 
to be {1, 2, • • • , n}. Then we have the decomposition 

1 
2 



E%" + Eg" = $[{£% + m ® (E?j, + E^,) - V=1(E% - E\t) % >/=I(i^, - £$*)], (8) 



-\{E, 



mn 
ij 



Eji ) 



1 



[m + e, 



•m \ 
Ik) 



l(E2 jf - E^,) + y/=l(E% - E%) ® {E$ r + &)]. (9) 



Equivalently we can picture the above decomposition as follows. We first view E^ n ± E^ n as 



:)■' 



an m x m block matrix (P s t), where P s t = except (s, t) = (k, I) or (/, k), and P^i = Pf k = E^y. 



Then we have E^ n + E^ n = E% 



E? f + E, 



l.k * 



also given by Eq.(8). 

Example 2. For / G [0, 1] consider the Werner state [19] 



Eyy. A simple computation will show that it is 



/ 1=1 
I 3 



1+2/ 1-4/ 

6 6 

1-4/ 1+2/ 



3 / 



(10) 



Then 
P 



+ 



1-f 

3 
1-/ 

3 
1-4/ 

12 



#1 + 



! + 2 /^4 



6 



(-^22 + £33) + 



!" 4 /^4 



() 



(^23 + ^32) + 



I"/ 



£ 



11 



l + 2f 
£11 <8> £11 + — tHH^ii ® #22 + £22 ® ^11) 



1-f 



[(E 12 + E 2 i) (£21 + E 12 ) - i(E 12 - E 2 i) <g> i(E 2 i - E 12 )] + -—^E 22 ® .E; 
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Example 3 For non-negative a, b, c consider the following positive semi-definite matrix 



Then we have 
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E u ® E n ® E u + E 22 <8> E 22 <8> E 22 

+t5'i2 <8 (S12 <£> S12 — iA\2 <8> iAw) — 7^12 <8 (£12 <8 M12 — i-Ai2 (8) S\ 2 ) 

+aE n ® E u <8 £22 + &-E11 <8 E 22 <8> En + c£ u E 22 ® E 22 

+-E 22 <g) En <8) En + tE 22 <8> En <8 E 22 + -E 22 <8) E 22 (8) En, 



where S, 



K) 



Eij + Eji and Aij 



E,-o — E 
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4 Separability of multipartite states 

As we note in the previous section that any Hermitian operator A on a tensor product space can be 
decomposed into a sum of tensor products of Hermitian operators: A = YH=i B\ ® Bf <g> ... <8> Bf. 
However the factors B\ are generally not density matrices on Hj as they may not be positive 
operators. To answer the question of separability of A one needs to study when each factor is 
non-negative. 

Let m(A) and M(A) denote the smallest and the largest eigenvalues of a Hermitian matrix 
A. We can transform the decomposition into another one so that the smallest eigenvalues are 
nonnegative as follows: 



A 



r 



E 2 ® ... ® Ef 



Y^ (b\ - m{B 1 i )Id 1 + m(B})IdA ® • • • (8 (b? - m{B?)Id n + m(B?)Id n 
i=i 
t 

^2 E- 1 <8> E- 2 (g) ... <g> E- n + g(A)Jdi Jcfc ® ••• ® M>, (12) 



one m(Bj? 



where B? are positive semi-definite Hermitian matrices on Hj, and each summand has at least 
but not all (i.e. at least one factor is the identity Id\ on Hi). 
Note that q(A) depends on the decomposition. We define the separability indicator of A, 
S(A) = max(q(A)) to be the maximum value of q{A) among all possible decompositions such as 
(12). The following result is quoted from [18]. 

Proposition 2. Let A = ^[ Ej (£> d be a density matrix on space H\ (8 H 2 . Then A is 
separable iff the separability indicator S(A) > 0. Moreover S(A) satisfies the following relation 
S(A) < m(A). 



Theorem 2. Let A = YH=i &i ® Cj be a Hermitian operator on the space H\ <g> H2, then q(A) 
is given by 

q(A) = Y m{B i )m(C i )+ Y m^M^i) 

m(Bi)>0,m(Ci)>0 m(B;)<0 

+ Y M(Bi)m(Ci)- Y m(Bi)m(Ci), (13) 

m(C;)<0 m(Bi)<0,m(Ci)<0 

and bounded by 



q(A) > M(A) - Y [(M(Bi) - m(Bij) (m(C 4 ) - m(Q)) + Mfm^Bi 

- m(m(Ci)B^ + M^miBAC^j - m{m{Bi)C^\. (14) 

Proof. For any Hermitian matrix P we define the operation P' by shifting with the minimum 
eigenvalue: P' = P — m(P)I. We can rewrite Eq. (12) 



A = ^ J B l , ®^+(m(C l ) J B l , -m((m(C J ) J B l , )/ m )®/„ 

+ I m ® (miBACi - m(m(Bi)Cl)I n ^ + q(A)I m ® J n , 

where 

?(A) = Y {rn(m(Ci)B^ + m(m(Bi)cfy + m{B t )m{Gi)) . 

i 

We observe that for any real s 

m(sP) = S+ ^ m{P) + S ~ ^ M(P), M(sP) = S + ^ M(P) + S ~ ^ m(P). (15) 
It then follows that 

q(A) = Y (m(m(Ci)Bfy + m(m(Bi)cfy + m(Bi)m(Ci 

i 

Y m{B i )m{C i )+ Y m{Bi)M{Ci)+ Y M(B t )m(CA 

m(Bi)>0,m(C z )>0 m(B l )<0,m(C* l )>0 m(Bi)>0,m(Ci)<0 

+ Y (m(Ci)M(Bi) + miBAMid) - m(Bi)m(C^ 

m(B l )<0,m(C i )<0 

Y m{B i )m{C i )+ Y rn(Bi)M(Ci) 

m(Bi)>0,m(Ci)>0 m(Bi)<0 

+ Y M(Bi)m(Ci)- Y rn(Bi)m(CA. 

m(Ci)<0 m(Bi)<0,m(Ci)<0 

Now we notice that for any matrix P and any real number r, M(P — rl) = M(P) — r, 
m(P — rl) = m(P) — r, from which it follows that 

M{B'A = M(B t - m(Bi)I m ) = M{BA - m(BA, 

MfacA&i - m{m{Ci)&AlJ) = M(m(Ci)B^j - m(m{Ci)B^ . (16) 



On the other hand it is well-known that M(A + B) < M(A) + M(B) (see [20]). Thus taking 
the maximum eigenvalues, we get 

r 

q(A) > M{A) - Y, [(M(Bi) - m(^)) (m(G) - m(Q)) + M(m{C l )B l 

i=\ 

- m(m(Ci)B^ + Af (m(Bi)Ci) - m(m{Bi)G, 

which completes the proof. 

In the last part of the proof if we take minimum eigenvalues we will get the known inequality 
m(A) > q(A) (using m(A + B) > m(A) + m(B)). 

We remark that the above lower bound is different from that in [18]. To better understand 
our lower bounds, we consider the special case when all factors are non-negative matrices, then 
m{m{Bi)Ci) = m(Bi)m(Ci) etc. Then it follows that 



q(A) > M(A) - Y [M(B t )M(C t ) - miBiJmid 
While the other extreme case is when all factors are negative, then 

r 

q(A) > M(A) - Y, [(M(Bi) - 2m(5«)) ( M (C*) " 2m(Q)) - m{B l )m{C l 
i=i 

When the factors Bi or G are not all non-negative, we have 
q(A) > M(A)- J2 [M{B i )M{C l )-m{B i )m{C i ) 

m(Bi)>0,rn(Ci)>0 

Y [(M{Bi) - 2m(Bi)) M(d) + m{B i )m{G 

m(Bi)<0,m(Ci)>0 

£ 

m(Bi)>0,m(Ci)<0 



(17) 



MiB^MiG) - 2m(Q)J +m{B i )m{C i 
M{Bi) - 2m(Bi)) (M{G) - 2m(C i )) - m(Bi)m(Ci) 



m(Bi)<nMCi)<n 
The above result can be generalized to multipartite states. 



(18) 



Theorem 3. Let A = Y^l Bi®Ci® Di be a density matrix on space H± <8> H2 <8> H3. Then A is 
separable if and only if the separability indicator S(A) > 0. Moreover S(A) satisfies the following 
relations 

S(A) < m(A). (19) 



r 
q{A) > M(A)-J2[( M ( B i)-m(B i )^(M(C i )-m(C i ))(M(D i )-m(D i 



i=i 



+ M^m(Bi)m(A)CiJ - m(m(B;)™(A)a) + M(ro(Ci)m(A)B. 

- m(m{Ci)m{Pi)B^ + M(m(B t )m{C z )D t ) - m(m{B i )m{C i )D l 



+ M(m(m(A)A)A - m(AMA)AJ - m(m(m(A)A)A - m(AMA)A 

+ M(m(m(A)A)A - m(A)m(A)A) - m[m{m{C i )B i )D l - m(C;)m(A)A 

+ M(m(m(A)A)Q - m(AMA)Q) - m(m(m(A)A)A - m(A)m(A)C, 

+ (M(m(A)A) - m(m(A)A)) (m(Q) - m(Q)) 

+ (M(m(C)A) - m(m(A)A)) (M(A) - m(A)) 

+ (M(m(A)Ci)-m(m(A)Ci))(M(A)-m(A))]- (20) 

The idea of the proof will be similar to that of Theorem 2 and is included in the Appendix. More 
generally we can use the same idea to give similar results for multi-partite cases. 

Theorem 4. Let A be a fc-partite mixed state on space H\ <S> A ® ■•• ® At, then A has a 
tensor decomposition into Hermitian operators in the form A = Y2i=i &1 ® Bf ® ••• ® A an d is 
separable if and only if the separability indicator S(A) > 0. Moreover S(A) satisfies the following 
relation 

S(A) < m{A). (21) 

When all factors are non-negative, we have 

r 

q (A) = J>(AV(A 2 )---™(A fc ) 

i 

r 

> M(A) - J2 [M(B})M(B?) ■ ■ ■ M(B?) - m(Bl)m(Bf) • • • m(B\ 



%=\ 



(22) 



5 Conclusion 

We have developed a criterion to judge whether a multi-partite density operator is separable. Our 
idea is first to decompose the density operator into a sum of tensor product of Hermitian operators. 
We give a new and practical way to decompose any Hermitian operator into tensor product of 
Hermitian operators in multi-partite cases. Unlike the numerical method [17] and the method 
of singular value decomposition [18] our new method is completely elementary and algebraic. 
Using the decomposition we can rewrite it into a tensor product of positive operators plus a 
scaler operator, which is called the separability indicator. The separability indicator provides a 
new mechanism to measure the quantum entanglement of the density operator. We derive some 
bound to estimate the scalar or separability indicator. Our inequalities are expressed in terms 
of eigenvalues of the summands, and in some case they are sufficient to tell if the separability 
indicator is non-negative, thus shows that the density operator is separable. As our method relies 
on how the operator is decomposed, it is usually difficult to compute the separability indicator 
exactly. We hope our estimates will shed more light on the separability problem. 
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6 Appendix 

Proof of Theorem 3. The idea of the proof is similar to that of Theorem 2. Recall the meaning 
of operation P' = P — m(P)I, and we have 



A = Y, B i® C i® D i = Y, B 'i® C i® D i + 12 { m ( D i) B i ~ m(m(Di)B^I m j ® CJ (8) /* 

j=l j=l 1=1 

r 

+ ^/m® (m(m(D t )B' i )C' l - m(ro(ro(A)B{)Cj)I„) ® 4 



i=i 



i=i 

r 

i=i 

r 

i=i 

r 



h 



i=\ 



® I n ® ifc 



i=i 

r 

+ J^ I m ® I n ® (m(m(C 4 )£-)^ - m(m(m(C i )Bl)D' l )I^ 
i=l 

r 

+ ^I m ®/„® [m^m^^-wf^^M^)^) 1 * + g(^)/ m ® J n ® J fc , (23) 



i=l 



where 



i=l i=l i=l 

r r r 

+ ^m(m(Ci)m(A)^) +^m(m(Si)m(A)C-) +^m(m(S i )m(C i )£>; 

i=l i=l i=l 

r 

+ £m(Bi)m(CiMA) 

i=l 
r r 

= ^m(m(m(A)^)Ci - m(A)m(Bi)Ci) +5^m(m(m(B i )C i )A - m(£;)m(Ci)A 
i=i i=i 

r r 

+ ^m(m(m(C7i)5i)A - m( J B i )m(C i )A) - ^m(m(A)^)m(a) 
»=i i=i 

r r r 

- Y,rn(m(B i )C l )m(D l )-Y,rn(HC i )B l )m(D i ) + Y,rn(B i )m(C l )m(D l ), (24) 



i=i 



i=i 



i=i 



where we have used similar identities like Eq. (16). Now we would like to consider eight possible 
signs of m(.Bj),m(Cj),m(A)j an d we use +,— ,+ to denote the subset {i\m(Bi) > 0, m(Ci) < 
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0,m(A) > 0} etc. to simplify the notation. By Eq. (15) it follows that 

q(A) = Y, ^(^)m(Q)m(A) 

+ Yl m ^) [M{C l )M{D i ) - M(Ci)m(Di) - m{C l )M{D l ) - m{C i )m{D l 

+ Yl m ( C ^ [M(Bi)M(Di) - m(Bi)M(Di) - M(^)m(A) - m{B i )m{D l 

+ Y m ^) [M(Bi)M(Ci) - m(Bi)M(Ci) - M(A)m(Q) - m(Bi)m(Ci) 

+ Y [m{Bi){M{C l )M{D l )-M{C i )m{D i )-m{C l )M{D l 

+ m {Ci) (M(Bi)M{Di) - m(^)M(A) - M(^)m(A)) 

+ Y [m(B i )[M(C l )M(D i )-M(C i )m(D i )-m(C l )M(D l 

+ m(A) (m(^)M(Q) - m^M^) - M(^)m(a)) 

+ ^ [m(C i )(M(^)M(A)-m(^)M(A)-M(^)m(A 

+ m(A) (M{Bi)M{Ci) - miB^Mid) - M{B l )m{C l 

+ Y (m(A)M(C i )M( J B i ) + m(^)M(Q)M(A) + m(Ci)M(A)M(B, 



2m(Bi)m(Ci)M(Di) - 2m{B i )m{D i )M{C i ) - 2m{C i )m(D i )M(B i 
+m{B i )m{C i )m(D i )). 



The above expression leads to an easy proof of the criterion: if A is separable, then all 
the factors are non-negative and S{A) > q{A) = s ^ ji m{Bi)m{Ci)m{Di) > 0. The converse is 
immediate. 

If we take minimum eigenvalues to the decomposition (12), we will get 



m(A) > Y m(B?)m(B?) • • • m(S>) + q(A) = q(A). 

i 

Next using similar identities as Eq. (16) we get identities like 
M(m(A)#- - m(m(A)B0/m)) = Af(ro(A)£i) - m(m(A)£i), 

Mfm(m(A)B t -)C< " m(m(m(A)£i)C-Kn 



(25) 



M(m(m(A)A)Q - m(Di)m{Bi)C, 
m(m(m{Di)Bi)Ci - m(Di)m(Bi)Ci 



Thus we have 



r 

M(A) < ^[(M(A)-m( J B,))(M(Q)-m(C l ))(M(A)-m(A 



i=i 



12 



+ M[ 

— m\ 

+ M 

+ M 

+ M 



(m(AMA)A) -m(m(5i)m(A)Ci) +M(m(Q)m(A)Bi 
(m(Ci)m(A)A) +M(m(A)m(Ci)A) -m(ro(A)m(Ci)A 
(m(m(A)Q)A - m(A)m(Ci)A) - m(m(m(A)Ci)A - ro(A)ro(Ci)A 
(m(m(Q)A)A - m(A)m(A)A) - m(m(m(A)A)A - m(A)m(A)A 
(m(m(A)A)C< - m(A)m(A)Ci) - m(m(m(A)A)C< - m(A)m(A)C, 
+ (M(m(A)A) - ro(m(A) A)) (^(Q) - m(C l )) 
+ (M (m(COA) - m(m(Q)A)) (m(A) - m(A)) 
+ (M(m(^)Q) - m(m(^)C 4 )) (m(A) - m(A))] + ?(A), (26) 

which completes the proof of Theorem 3. 

We remark that when all factors are non-negative matrices, then it follows that 



(27) 



q(A) > M(A) - ^ \M(Bi)M(Ci)M(Di) - m(A)m(C«)m(A 

j=l 

While the other extreme case is when all factors are negative, then 

r 

q(A) > M(A)-J2 [(Af(A)-3m(A))(Af(C i )-3m(A))(M(A)-3m(A))-m(A)m(A)m(£>. 



i=i 



Proof of Theorem 4. The proof is by an easy induction as those of Theorems 2 and 3. 
Some details are already offered in Eqs (25) and (26). 
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